Entropic force of polymers on a cone tip 
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We consider polymers attached to the tip of a cone, and the resulting force due to entropy loss on 
approaching a plate (or another cone). At separations shorter than the polymer radius of gyration 
R g , the only relevant length scale is the tip-plate (or tip-tip) separation h, and the entropic force 
is given by F = AksT '/h. The universal amplitude A can be related to (geometry dependent) 
correlation exponents of long polymers. We compute A for phantom polymers, and for self-avoiding 
(including star) polymers by e-expansion, as well as by numerical simulations in 3 dimensions. 
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Single molecule manipulation [J-Q using techniques 
such as atomic force microscopy (AFM) Q, micronee- 
dles 0] , optical [1, H| and magnetic [T3| tweezers enable 
extremely detailed study of geometry and forces in long 
polymers. The positional accuracy of AFM tip JEHU can 
be as good as few nm, while the forces of order of 1 pN 
can be measured, and measurements can be carried out 
in almost biological conditions [HI, [l]| . These enhanced 
sensitivities bring us to the range where entropic forces of 
long polymers in a solvent can be significant even when 
the deformation of the polymer is relatively slight. While 
the main thrust of the experimental research is extrac- 
tion of specific information from the force-displacement 
behaviors of the polymers, certain features are indepen- 
dent of the microscopic details [l4j . but depend on the 
probe shape, as discussed in this work. 

We consider an idealized set-up in which a polymer is 
attached to the tip of a solid cone. The cone approaching 
a plate (or another cone) exemplifies a geometry in which 
the only (non-microscopic) length scale is provided by 
the tip-plate (or tip-tip) separation h. Fluctuating poly- 
mers typify self-similar variations at scales intermediate 
between microscopic (persistence length a) and macro- 
scopic. The latter is set by the radius of gyration which 
grows with the number of monomers through the scal- 
ing relation R g oc N v . Thus when a cone-tip-attached 
polymer approaches a plate, at separations a -C h <C R g 
the only relevant length scale is h, and on dimensional 
grounds, the force due to loss of entropy must behave as 



F = A 



(1) 



Such a force law should apply to all circumstances where 
the separation provides the only relevant length scale. 
The amplitude A will depend on geometric factors such 
as the opening angle of the cone O (and if tilted, on the 
corresponding angle). One could presume, that the di- 
mensionless amplitude may also depend on microscopic 



properties on the polymer. However, in case of cone-tip- 
polymers we shall demonstrate that the amplitude A can 
be related to universal (and shape dependent) polymer 
exponents. The simple force law of Eq. (JXJ) follows easily 
from various polymer scaling forms (see, e.g. the deriva- 
tion below) such as in Refs. [15l - [r7| . and should be part of 
polymer lore. Surprisingly, we could not find an explicit 
reference to it in any of the standard polymer textbooks. 

A polymer attached to the tip of an AFM is approx- 
imated as linked to the apex of a cone as depicted in 
Fig. [TJ,. With the cone far away from a plate (h ^> R g ), 
the number of configurations of the polymer grows with 
the number of monomers as 



JV C = bz N N lc( - e ^ 



(2) 



where the effective coordination number z, as well as the 
pre-factor 6, depend on the microscopic details, while the 
'universal' exponent j c only depends on the cone angle. 
When the cone touches the plate as in Fig. [Dp, the num- 
ber of configurations is reduced to A/" cp with the same 
form as Eq. @, but with a different exponent 7 cp (0). 
We shall henceforth use the exponent subscript 's' (as in 
7 S ) to refer to the above cases, with "s=c" for cone and 
"s=cp" for cone+plate; the absence of a subscript (as in 
7) will signify a free polymer. The work done against the 
entropic force in bringing in the tip from afar to contact 
the plate can now be computed from Eq. |T]) as 



W 



"dhA^^ =Ak B T\n^- = Avk B T\nN. (3) 
h a 



The work can also be computed from the change in free 
energies between the final and initial states, due to the 
change in entropy, as 



AT 



-TAS = TS C - TS C 



k B T( lc - lcp )lnN, (4) 

where the entropy S = — fcslnA/" was computed from 
Eq. |2]). By equating W and AT we find 



A 



7c - 7c 



(5) 
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the final result obtained from the scaling law 7 S = 
(2 — ris)^, where 77 characterizes the anomalous decay of 
correlations (~ \/r d ~ 2+r >). 
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FIG. 1. (a) Polymer attached to the tip of a solid cone with 
apex semi-angle <d (configuration "c"); positions are described 
by the spherical coordinates r, 9 and azimuthal angle cj> (not 
shown), (b) The tip is at a distance h <C R g from the plate, 
(c) The tip touching the plate (configuration "cp"). 



In the above discussion we have assumed that the only 
interaction between the polymer and the surfaces is due 
to hard-core exclusion. Attractive interactions between 
the polymer and surface will introduce temperature de- 
pendent corrections, and an additional size scale. Weak 
interactions are asymptotically irrelevant, but strong in- 
teractions may lead to a phase in which the polymer is 
absorbed to the surface, where the entropic considera- 
tions presented here will no longer be appropriate. 

We have derived Eq. Q for the specific case of a cone 
and a plate. However, this equation can be applied to 
any situation where in the limits of large and vanishing 
separations we have scale-free shapes. In particular, the 
same rule can be applied when a linear or star polymer is 
brought from infinity to a contact with a repulsive plane, 
since both extremes have no length-scale, leading to exact 
expressions for the force constant (l6| . However, if the 
AFM tip is slightly rounded, an additional length scale 
is introduced, and one may expect a non-trivial crossover 
between various regimes jl8j ]. 

We have thus reduced the computation of the force 
to calculation of correlation functions. In the absence 
of self-avoidance and other interactions, correlations of 
the so-called ideal or phantom polymer (henceforth de- 
noted by subscript 0) are the same as a free-field theory, 
satisfying (at scales shorter than R g ) the Laplace equa- 
tion V 2 G (r,r') = -<5(r - r') With one point at 
a short distance a from the cone, correlations behave as 
G (a,r) ~ a 710 /r d - 2+r >°ty(9), where r and 6 denote the 
distance of r from the tip, and the angle of r to the cone 
axis. The change in scaling from a free phantom poly- 
mer is captured by (a/r) Vo , and ^f(9) is a dimensionless 
function depending only on the polar angle due to the 
symmetry of the geometry. (We consider a generalized 
cone in d spatial dimensions characterized by a single 
polar angle 9, and d — 2 azimuthal angles tf>, ip, ■ ■ ■ .) Sub- 
stituting the above form in the Laplace equation, we find 
that the exponent t]q satisfies 
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77o)*(0) = O, 
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FIG. 2. (Color online) The exponent ry s o for ideal polymers 
in d =2 (dot-dashed), 3 (solid), 4 (dashed) for cone ("s=c") 
of angle B (bottom curves), and "s=cp" (top curves). 



with an appropriate boundary condition on "J. For an 
isolated cone, the function ^ must be positive and regu- 
lar outside the cone, with d^/d9\$ =7I — to avoid a cusp 
on the symmetry axis, and ^(0) = on the cone surface. 
For the cone+plate, the appropriate solution is positive 
and vanishes both at 9 = and 9 = ir/2. The first case 
was considered by Ben-Naim and Krapivsky [l9[ in con- 
nection with diffusion near an absorbing boundary [20| . 
and we follow these derivations. The solution in general 
d requires the use of associated Legendre functions, but 
simplifies in a few cases described below. 
• For d = 2, Eq. (JBJ) reduces to + r) 2l $> = 0; solved 
by linear combinations of sm(r] 9) and cos(t]q9) to yield 



2(tt-0)' 



and 



2tt 



VcpQ 



20 



(7) 



Both results (depicted in Fig. [5J go to a finite value as 
— > 0, reflecting the strong reduction in configurations 
due to the remnant (barrier) line, and rjo — > oo when the 
boundaries confine the polymer to a vanishing sector. 
• For d — A, the substitution ^ = u/ sm.9 simplifies 
Eq. ([6]) to u" + (rjo + l) 2 u = 0, solved by a linear combi- 
nation of sm[(j7o + 1)9] and cos[(??o + 1) 







7T-0' 



and rj, 



], and we find 
f 20 



cpO 



7T-20 



(8) 



depicted by the bottom and top dashed lines in Fig. [2] 
The cone exponent rj c o vanishes linearly with — a needle 
in four dimensions is invisible. 

• For d — 3, Eq. (jSJ) is solved by a linear combination 
of regular Legendre functions P rj0 (cos9) and Q rj0 (cos9). 
The resulting exponents (which cannot be cast as simple 
functions), are plotted as solid lines in Fig. [2j The ex- 
ponent ?7co(0) vanishes with the angle, but both curves 
approach their limiting value as — ¥ with infinite slope 
via a logarithmic singularity (~ 1/| ln0|). 
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• For all d, the cone becomes a plate for = tt/2. Cor- 
relations with one point approaching a surface are easily 
obtained by the method of images [2l| leading to T] c o = 1, 
which is clearly seen in Fig. [5] For 3 < d < 4, both expo- 
nents approach their limiting value when — > as Po 
with po = d — 3: 



VcpO 



= 1 



r(l-e/2) 
0Fr(l/2-e/2) 
4r(2 - e/2) 



1 



and 



/7pr(l/2- e/2) 



1 



(9) 



where d — A—e. These equations reflect the fact that the 
two dimensional phantom polymer will not intersect the 
remnant line in d > 3. For d < 3, the limiting value is 
different from the case without any cone indicating the 
finite probability of intersection of the polymer with the 
line. 

Universal aspects of swollen (coil) polymers with short- 
range interactions can be modeled by a self-avoiding walk 
(SAW). In d = 3, a SAW in empty space has exponent 
7 « 1.158 2l. while in d = 2, 7 = 43/32 [H (for ideal 
polymers 7 = 1 at any d) . There are a number of results 
regarding 7 S for polymers confined by wedges in 2D and 
3D [24- 29]: A SAW anchored at the origin and confined 
to a solid wedge (in 3D) or a planar wedge (in 2D) has an 
angle-dependent 7 S that diverges as the confining angle 
vanishes. Studies of a polymer attached to the tip of a 
2D sector in 3D, and to the apex of a cone have been per- 
formed [3Cj ■ Extensive analytical [U [H| and numerical 
[H, [2t| studies of SAWs anchored to a solid plate in 3D 
find 7 S = 71 « 0.70 H|, or 0.68 [33[. For ideal polymers 
71 = 1/2 for any d. We are unaware of specific results 
for the geometries depicted in Fig. [TJ and report below 
our numerical and analytical estimates. 

We performed numerical simulations for SAWs of 
lengths N = 16, 32, ... , 1024 on a cubic lattice. The 
SAWs were generated by a dimerization method [U, |35[ 
in which an unbiased iV-step SAW is created by attempt- 
ing to join two 7V/2-step SAWs previously obtained by 
the same algorithm. We generated 10 8 SAWs for each 
N, each of which is attached to the origin and checked 
whether it touches a cone (or cone+plate). The probabil- 
ity that an JV-step SAW does not intersect the excluded 
space is the ratio of permitted number of walks to the to- 
tal number of SAWs, i.e. p N = N s /M ~ N^/N^- 1 = 
N AlB . The ratio pn/p2N — 2 A7b is then used to esti- 
mate A7S for each sequential pair of ./Vs. The results 
are extrapolated to N —> 00 by plotting the estimates 
versus l/y/N; errors are due to both finite sample size 
and insufficiently large N. The estimated exponents (full 
symbols in Fig. [3]) are rather close to the results for ideal 
polymers. This apparent proximity is somewhat mislead- 
ing as A7S = 7 — 7s masks an appreciable shift in both 
7 and 7s- For example, for a SAW near an isolated cone 
we expect A7(0 = tt/2) = 7 — 72 =0.44-0.48, very close 
to ideal polymer result of 1/2; our numerical estimate at 
= tt/2 is 0.477 ±0.004. 




FIG. 3. (Color online) Dependence of the exponent difference 
A7 S = 7 — 7 S on angle Q for the cone+plate (top curves), and 
an isolated cone (bottom curves). The solid curves are the 
results for an ideal polymer, while the data points represent 
numerical results for SAWs in the same geometry. All error 
bars are estimates of the N — > 00 extrapolation error. The 
dashed line depicts prediction of the e-expansion in Ref. [30l ]. 



The dashed line in Fig. [3] depicts the result [30j of an 
e = (4 — d) expansion in which the 2-dimensional surface 
of the cone is treated as a weakly repulsive potential. 
The lowest order result, A7 C = (3e/8)sin0, captures 
some features of the numerics for < tt/2, but fails dis- 
mally for > 7r/2 (since the polymer simply jumps to 
the larger space inside the cone). It also fails to cap- 
ture the correct behavior as — > 0. A better approach 
is to exclude the interior of the cone, and to this end 
we follow the work of Cardy [25j, [26J for the wedge ge- 
ometry. The corresponding computations for a cone are 
more complicated, and we relied upon recent results on 
the electrodynamic Casimir force on a cone [3||. 

For an e-expansion we need the full non-interacting 
Green's function in 4-dimensions, going beyond the large 
separation asymptotic form obtained via Eq. (For a 
general e-expansion, one has to know the Green's func- 
tion in 4 — e dimensions, however, we carry out this ex- 
pansion only to first order in e. Since the strength of the 
interaction is linear in e, it suffices to obtain the Green's 
function in 4-dimensions.) Following Ref. [3o| . this is 
given by 



G (x,x') 



£ (-i)'tt r( Pfc +z+i) 



-1-1/2/ 
' gfc-Vg * 



Vsin 8' 



(10) 



where the sum is over the triplet of integers k > 0, I > 0, 
and — I < m < +1. The important exponent pk is the k-th 



root of the transcendental equation P 



-1-1/2 



(-COS0) 



0. This Green's function is broken up in radii (r < and 



4 



r>), as appropriate to a polymer in the presence of a 
cone with one endpoint close to the tip and the other 
end far away. 

The interacting Green's function at first order is ob- 
tained by subtracting polymer configurations that self- 
intersect, forming an intermediate loop, as 

G x = G - u[ d 4 x"G (x,x")G r (x",x")G a (x",x'). (11) 



In the language of quantum field theory, the first term is 
the "free" Green's function, albeit in the presence of ex- 
ternal boundary conditions, while the second term is the 
one-loop correction. Note that the intermediate Green's 
function is regularized by subtracting the Green's func- 
tion in empty space, hence the superscript r. In the above 
equation, u is the strength of the self-avoiding interac- 
tion, which is ultimately set to its fixed point value of 
u* = 2n 2 e [111, [IEEm!- To calculate the scaling proper- 
ties for r 3> a, it is sufficient to include only the first term 
(k = 1, 1 = m = 0) from Eq. (|T0|) in the non-loop prop- 
agators (as in the wedge computation by Cardy [25|, [26j , 
higher-order terms give corrections in higher powers of 
a/r). The intermediate loop, however, can be of any 
size requiring the entire sum (albeit regularized by sub- 
tracting the result for 9 = to remove an unrelated 
divergence). The loop correction is an integral over the 
whole space. The integration over angles 4>" and ip" is 
trivial due to symmetry and yields 1(1 + 1) when summed 
over all spherical harmonic functions of degree I; in the 
integration over the radius r" we seek a logarithmic con- 
tribution which is exponentiated in the end. Similar to 
Cardy's analysis, such a logarithm comes only from the 
region r < r" < r' . There remains an integral over the 
polar angle 9", but this is cumbersome as we have to 
find the roots to the transcendental equation noted be- 
fore. We shall not dwell further on the e-expansion for 
a general opening angle. Instead, we focus on the limit 
of a sharp cone as 9 — > 0, where the leading singularity 
comes only from the I = channel of the loop Green's 
function. This is because a sharp cone couples only to the 
lowest spherical partial wave. A careful integral over 9" 
and summation over all k's gives the radial dependence 
of the Green's function as 



G x cx 



r »icO 



r 1 
1 + eln- — 91n9 
r' V 4-7T 



.166 



(12) 



where the term proportional to e in the bracket is the 
one-loop correction to the Green's function. In d = 4 — e 
dimensions, r) cQ is given by Eq. Now we can expo- 
nentiate the radial logarithm to obtain the renormalized 
exponent, 



T)c = VcO + 



eine 

47T 



.16 9 



(13) 



Note that the loop-correction vanishes logarithmically 
with the angle, as there is no first order in e contribution 



to r\ in empty space. Repeating the same procedures for 
the cone-plate we find, as 6 — » 0, 



Vcp = VcpO + 



1 eine 



.66 6 



(14) 



where the exponent 77 cp o is given by Eq. (j^J). In this case, 
the loop-correction goes to — e/8 for 9 = due to the 
presence of the plate. Equations (|T5j) and (|14p suggest 
interpreting the logarithmic corrections as signatures of 
a power-law for 9 — > 0. Expanding rj c0 and r/ cp0 to the 
first order in e, we find another In 9, originating from the 
expansion in 4 — e dimensions of the phantom polymer 
(as opposed to the perturbative terms from the one-loop 
computation). Summing both contributions, we obtain 



9 / 3 

r) c = — 1 e In 9 — .06e I , and 

7T \ 4 



»fcp = 1 - I + ^ (l - ^eln6 - .86e 



(15) 



We can then recast the approach of the exponents to 
their limiting values as a power-law 9 P with p = 1 — 3e/4 
in place of po = d — 3 = 1 — e for the phantom poly- 
mer (Fig. [2J. We may interpret this result as follows: 
self-avoidance swells the polymer, reducing its fractal di- 
mension from 2 to v^ 1 = 2 — e/4 to lowest order. The 
dimensionality of the intersection of such an object with 
the remnant line is <i — 1 — v^ 1 = 1 — 3e/A = p. We leave 
this observation as a conjecture for future studies. 

From a practical point of view, the amplitude A = 
Vcp ~ Vc is typically a number of order unity. Thus the 
force in Eq. (fTJ at room temperature is roughly around 
0.1 pN at 0.1 /im separation; at the margins of possi- 
ble measurements by current force apparatus. The force 
can be increased by attaching several polymers to the 
tip. For ideal (phantom) polymers the force is enhanced 
by /, the number of polymers, while interactions mod- 
ify this conclusion. The interactions are of two kinds: 
self-interaction of a single arm (intra-arm) and interac- 
tions between two different arms (inter-arm) . The former 
contribution is computed above to first order in e and 
should be multiplied by / for an /-arm polymer. The 
latter, however, gives rise to a new term multiplied by 
/(/ — l)/2, the number of pairs. Thus to first order in e, 

i] = fr)o + fr]i + ^2 — r l e wnere ^0 denotes the exponent 
of a single phantom polymer, and r\i and r\ e are those 
of intra-arm and inter-arm interactions respectively. In 
the absence of a cone, the set-up is similar to the widely 
studied case of star polymers [37| . We carried the corre- 
sponding e-expansion with cone and cone+plate. Here, 
we just quote the final result after appropriate exponen- 
tiations. We find (in the limit of sharp cones), a force 
amplitude 



Af) 



e 

X -8 



3 



.80 



£0 



0l-3e/4 _ 

(16) 
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Interactions amongst the polymers thus reduce the force 
amplitude per polymer. 

Note that the same exponent p = 1 — 3e/4 governs the 
approach to a finite limit in Eq.(fl6|) as — > 0. Another 
experimental set-up is a cone-tip-attached polymer ap- 
proaching another cone. The entropic force is of course 
much smaller in this case, and we have verified that it 
also vanishes with cone angle as O p , bolstering the con- 
jecture that this is a universal exponent related to the 
needle geometry. 

In summary, we propose that polymers exert an en- 
tropic force AksT/h on a cone tip, with a 'universal' am- 
plitude A dependent on geometry, interactions, and num- 
ber of polymers. We conjecture that the singular form 



of the amplitude on vanishing cone angle is described by 
a new exponent. There are many set-ups where a simi- 
lar force law is expected on the basis of scaling at length 
scales shorter than an appropriate correlation length. 
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